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Abstract 

In models with an extra U(l)' gauge boson family non-universal couplings to the weak eigen- 
states of the standard model fermions generally induce flavor-changing neutral currents. This 
phenomenon leads to interesting results in various B meson decays, for which recent data indicate 
hints of new physics involving significant contributions from b —* s transitions. We analyze the B s 
system, emphasizing the effects of a Z' on the mass difference and CP asymmetries. 
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I. INTRODUCTION 



The study of B physics and the associated CP violating observables has been suggested as 
a good means to extract information on new physics at low energy scales 
Since B-B mixing is a loop- mediated process within the standard model (SM), it offers 
an opportunity to see the footprints of physics beyond the SM. The currently observed 
AM d = 0.489 ± 0.008 ps _1 M and its mixing phase sin 23 = 0.736 ± 0.049 extracted from 

Q n 

the Bd — > J/ipKs mode jj9j agree well with constraints obtained from other experiments [10]. 



■ exper 



However, no such information other than a lower bound AM S > 14.5 ps [11] is available 
for the B s meson yet. 

Based upon SM predictions, AM Bs is expected to be about 18 ps" 1 and its mixing phase 
<f) s only a couple of degrees. In contrast to the Bd system, the more than 25 times larger 
oscillation frequency and a factor of four lower hadronization rate from b quarks pose the 
primary challenges in the study of B s oscillation and CP asymmetries. Since the B s — > J/i/j(f> 
decay is dominated by a Cabibbo-Kobayashi-Maskawa (CKM) favored tree-level process, 
b — > ccs, that does not involve a different weak phase in the SM, its asymmetry provides 
the most reliable information about the mixing phase cj) s . 

Although new physics contributions may not compete with the SM processes in most of 
the b — > ccs decays, they can play an important role in B s -B s mixing because of its loop 
nature in the SM [12j. In particular, the mixing can be significantly modified in models in 
which a tree-level b — > s transition is present. Thus, measurement of the properties of B s 
meson mixing is of high interest in future B physics studies as a means to reveal new physics 
13, [l^. Since the current B factories do not run at the T(5S) resonance to produce B s 



mesons, it is one of the primary objectives of hadronic colliders to study B s oscillation and 
decay in the coming years 

HQ. 

Flavor changing neutral currents (FCNC) coupled to an extra U(l)' gauge boson arise 
when the Z' couplings to physical fermion eigenstates are non-diagonal. One way for this to 
happen is by the introduction of exotic fermions with different U(l)' charges that mix with 



the SM fermions 



18. 
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211 ] as occurs in E 6 models. In the E e case, mixing of the 



right-handed ordinary and exotic quarks, all SU{2)l singlets, induces FCNC mediated by 
a heavy Z' or by (small) Z-Z' mixing, so the quark mixing can be large. Mixing between 
ordinary (doublet) and exotic (singlet) left-handed quarks induces FCNC mediated by the 



2 



SM Z boson 2l|. We will also allow for this possibility, but in this case the quark mixing 
must be very small. 

Another possibility involves family non-universal couplings. It is well-known that string 



models naturally give extra U( . 
couplings to the SM fermions (24 , 



grou ps, at least one of which has family non- universal 
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bzi- 



Generically, the physical and gauge eigenstates 
do not coincide. Here, unlike the above-mentioned Eq case, off-diagonal couplings of fermions 
to the Z' boson can be obtained without mixing with additional fermion states. In these 
types of models, both left-handed and right-handed fermions can have family non-diagonal 
couplings with the Z', while couplings to the Z are family diagonal (up to small effects from 
Z — Z' mixing). 

The Z' contributions to B s -B s mixing are related to those for hadronic, semileptonic, 
and leptonic B decays in specific models in which the diagonal Z' couplings to qq, £ + £~, 
etc. are known, but are independent in general (46 1. We have found that in specific models, 
B s -B s mixing effects can be significant while being consistent with the other constraints; 
these results will be presented elsewhere. In the present paper, we will treat the mixing in 
a model-independent way. 

between 
which 



Recently, we have studied the implications of a sizeable off-diagonal Z' coupling 



the bottom and strange quark in the indirect CP asymmetry of B — > <t> K<j de cay 
appears to show a significant deviation from the SM prediction [5j, |6j, |29|, [3]. Here we 



extend our analysis to B s -B s mixing where the Z' contributions also enter at the tree level. 
Applications to the B — > irK anomaly are under investigation j^jj. 

The paper is organized as follows. In Section |Hl we review the basic formalism of B s - 
B s mixing. In Section IIII1 we evaluate AM S in the SM. In Section IIV1 we include the 
Z' contributions, allowing both left-handed and right-handed couplings in the mixing, and 
study their effects on observables. Our main results are summarized in Section IVl 



II. B,-B H MIXING 



In the conventional decomposition of the heavy and light eigenstates 

\B s ) L =p\B s ) + q\B° s ) , 
\B s ) H =p\B° s )-q\B° s ) , 



(1) 
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the mixing factor 



Q 



has a phase 



Mf 2 M * 



<>f M = 2 arg( V t6 V t *) = -2\ 2 fj ~ -2° , sin 20f M ~ -0.07 , (3) 



where the theoretical expectation r^f ^ Mf 2 M is used. The approximate formula Eq. (J2J) 
receives a small correction once Tf^ 1 is included. Model independently, this only shifts 4> s is 
at the few percent level. With errors on A and fj included, we have the SM expectation that 
sin20f M ~ -0.07 ±0.01. 

The off-diagonal element of the decay matrix, rf^, is evaluated by considering decay 
channels that are common to both B s and B s mesons, and M\i is the off-diagonal element 
of the mass matrix. Due to the CKM enhancement, Tf^ 1 is dominated by the charm-quark 
contributions over the up-quark contribution in a box diagram. Unlike the Kaon system, 
rf^ is much smaller than Mf 2 M for B mesons because the former is related to the B meson 
decays and set by the scale of its mass, whereas the latter is proportional to rn? t . We can 
safely assume that Ti2 is not significantly modified by new physics because Ti2 receives major 
contributions from CKM favored b — *> ccs decays in the SM, and the SM result I\ 2 <C Myi 
is unlikely to change. 

The mass difference of the two physical states is 

AM S = M H — M L ~ 2|M 12 | . (4) 

The width difference is 

2Re(M 1 * 2 r 12 ) 

AT = T H -T L = — — = 2|r 12 | cos9 , (5) 

where the relative phase is 6 = arg(Mi 2 /ri 2 ). Since Ti 2 is dominated by the contributions 
from CKM favored b — > ccs decays, we have 6 = arg ((V t bV* s ) / (V c bV* s )) ~ n, and thus 
Ar ~ — 2 1 1 is negative in the SM. Although Ti 2 is unlikely to be affected by new physics, 
the width difference always increases as long as the weak phase of M 12 gets modified |32j. 
The observability of B s -B s oscillations is often indicated by the parameter 

AM, 



X u 



(6) 



4 



where T s = (4.51 ± 0.18) x 10 13 GeV, converted from the world average lifetime r s = 
1.461 ± 0.057 ps [8]. The expected large value of Xojs a challenge for experimental searches. 
Currently, the result from all ALEPH \3^, CDF [3J, DELPHI 35], OPAL |3JJ, and SLD 
37| studies of AM S with a combined 95% confidence level (CL) sensitivity on AM S of 18.3 
ps -1 gives jll | 

AM, > 14.5 ps" 1 , and x s > 20.8 . (7) 

It is also measured that m Ba = 5369.6 ± 2.4 MeV Q and AJ\/T S = -0.16^(|Ar s |/r s < 
0.54) (with the 95% CL upper bound given in parentheses [rjj) consistent with recent next- 
to-leading-order (NLO) QCD estimates [2^. In comparison, the B d system has m Bd = 
5279.4 ±0.5 MeV, AM d = (0.489 ± 0.008) ps" 1 , x d = 0.755 ±0.015, and r Bd = 1.542 ±0.076 
ps 



III. AM S IN THE SM 

The \AB\ = 2 and |AS"| = 2 operators relevant for our discussions are: 

LL = [s 1 ,(l- l5 )b}[sY(l-l,)b} , 
0^ R =[s 1 ,(l- l5 )b}[sY(l+l 5 )b} , 

Of* = [3(1 ~ 75)6] [5(1 + 75)6] , 

RR = [57 M (1 + 75)6] [«7"(1+ 75)6] • (8) 

Because of the V — A structure, only the operator LL contributes to B s -B s mixing in 
the SM. The other three operators appear in the Z' models because of the right-handed 
couplings and operator mixing through renormalization, as considered in the next section. 
In the SM the contributions to 

M!r ^ ^-(B° s \ntm) (9) 

ivn Bs 

are dominated by the top quark loop. The result, accurate to NLO in QCD, is given by 



Mf 2 M = I ^M 2 Vm B Ji s (^K t :) 2 r /2i3 5o(a ;t )[a s (m b )]- 6 / 23 
where x t = (m t (m t )/M w ) 2 and 



1 + ^T^ 



B LL {m b ) , (10) 



Ax — llx 2 + 3a; 3 In a; 

S0{X) = 4(1 -s)* W^? • ( } 



Using m t (m t ) = 170 ± 5 GeV, we find S (x t ) = 2.463. The NLO short- distance QCD 
corrections are encoded in the parameters r]2B — 0.551 and J5 ~ 1.627 [39]. The bag 
parameter B LL (p) is defined through the relation 



(B s \0 LL \B s ) = -mljlB LL ^) . 



(12) 



In the following numerical analysis, we will use Gp = 1.16639 x 10~ 5 GeV" 2 and M w 
80.423 ± 0.039 GeV [8], and write the SM part of AM S as 



AM! 



SM 



1.19 



0.04 



ft 



230 MeV 



B LL (m b ) 
0.872 



x 10" 11 GeV . 



(13) 



Current lattice calculations still show quite lar ge e rrors on the hadronic parameters 



230 ±30 MeV and B LL { 



m b ) = 0.872 ± 0.005 



0. 



4. 



421 ] . However, the ratio 



fB d \fB Bd 



(14) 



can be determined with a much smaller theoretical error, where Bg is the renomalizat ion- 
independent bag parameter for the B q meson (q = d,s). Therefore, the error on AM S within 
the SM can be evaluated by comparing with AM^, i.e., 

(1-A 2 ) 2 



AM S SM = AM| M £ 



SM t mB ° 



(15) 



m Bd X 2 [(l-p)2+jf] 

Using the measured values of the Wolfenstein parameters 43j A = 0.2265 ± 0.0024, A = 
0.801 ± 0.025, p = 0.189 ± 0.079, and r/ = 0.358 ± 0.044 [loj, ^ = 1.24 ± 0.07 |44j, and the 
mass parameters quoted above, we obtain the SM predictions 

AMf = (1.19 ±0.24) x 10" 11 GeV 
= 18.0 ±3.7 ps" 1 , 
26.3 ±5.5 . 



x 



SM 



(16) 



As noted above, the central value of x s is slightly larger than the current sensitivity based 
upon the world average. Recent LHC studies show that with one year of data, AM S can be 
explored up to 30 ps _1 (ATLAS), 26 ps _1 (CMS), and 48 ps _1 (LHCb) (corresponding to 
x s up to 46, 42, and 75); the LHCb result is based on exclusive hadronic decay modes [l6j. 
The sensitivity of both CDF and BTeV on x s can also reach up to 75 using the same modes 
[lfi| . for a luminosity of 2 fb^ 1 . The sensitivity on sin20 s is correlated with the value of x s , 
and it becomes worse as x s increases. A statistical error of a few times 10~ 2 can be reached 
at CMS and LHCb for moderate x s ~ 40 Q- 



IV. Z 1 CONTRIBUTIONS 



For simplicity, we assume that there is no mixing between the SM Z and the Z' (small 



mixing effects can be easily incorporated 17[). A purely left-handed off-diagonal Z' coupling 
to b and s quarks results in an effective \AB\ = 2, \AS\ = 2 Hamiltonian at the Mw scale 

HS = % S %Ae^ (mt) , (17) 

where #2 is the U(l)' gauge coupling, g% = e/(sin8w cosOw), Mz' is the mass of the Z' ' , and 
B^ b is the FCNC Z' coupling to the bottom and strange quarks. The parameters pi and the 
weak phase 0l in the Z' model are defined by the second equality. Generically, we expect 
that g-z/ Q\ ~ 1 if both U(l) groups have the same origin from some grand unified theory, and 
M z /M Z i ~ 0.1 for a TeV-scale Z' . If \B^ b \ ~ |V^VJ*|, then an order-of-magnitude estimate 
gives us pi ~ O(10~ 3 ), which is in the ballpark of giving significant contributions to the 
B s -B s mixing. The Z' does not contribute to Ti2 at tree level because the intermediate 
Z' cannot be on shell. After evolving from the My/ scale to m^, the effective Hamiltonian 
becomes 



niZ' 



Gp 
V2 



a s (m b ) - a s (M w ) 
1+ 4^ Js 



R^p 2 L e 2i ^O LL {m b ) , (18) 



where R = a s (Mw) / 'a s (m(,) . Although the above effective Hamiltonian is largely suppressed 
by the ratio (g2Mz)/(g\Mz'), it contains only one power of Gf in comparison with the 
corresponding quadratic dependence in the SM because the Z'-mediated process occurs at 
tree level. 

The full descrirjtion of the running of the Wilson coefficient from the M\y scale to rrib 
can be found in We only repeat the directly relevant steps here. The renormalization 
group equation for the Wilson coefficients C, 

C = 1 T (g)C(p), (19) 



d In p 

can be solved with the help of the U matrix 

C(p) = U(p,M w )C(M w ) , (20) 

in which 7 T (g) is the transpose of the anomalous dimension matrix 7 (g). With the help 
of dg/dliap = 13(g), U obeys the same equation as C(p). We expand j(g) to the first two 



terms in the perturbative expansion, 



Of. 



7-^ + 7 



(i) 



4tt/ 



(21) 



(22) 



To this order the evolution matrix U(p, m) is given by 

U(p,rn) = (l+ U^m) (l - *J 

where f/^ ^ is the evolution matrix in leading logarithmic approximation and the matrix J 
expresses the next-to-leading corrections. We have 

/ 



m) 



V 



v- 1 



(23) 



D 



where V diagonalizes 7 ( -°^ T , i.e., 7^ = V~ 1 7(°) T V, and 7^ is the vector containing the 
diagonal elements of the diagonal matrix 7^. In terms of G — V~ 1 'y( 1 " >T V and a matrix H 
whose elements are 



Ha 



(0) _A 

2/?o 2 2/5 + 7 



(°) __ J°) 
13 



(24) 



the matrix J is given by J = VHV^ 1 . 

The operators LL and Ri? do not mix with others under renormalization. Their Wilson 
coefficients follow exactly the same RGEs, where the above-mentioned matrices are all simple 
numbers. The factor 

1 R 6 / 23 (25) 



a s (m b ) - a s (M w ) 
1 H ; J 5 



4tt 

in Eq. (|T8|) reflects the RGE running. On the other hand, Of R and 0% R form a sector that 
is mixed under RG running. Although the Z' boson only induces the operator 0\ R at high 
energy scales, 0\ R is generated after evolution down to low energy scales and, in particular, 
its Wilson coefficient C% R is strongly enhanced by the RG effects 

With contributions from both the SM and the Z' boson with only left-handed FCNC 
couplings included, the B s mass difference is 



AM, 



AM - M f 1 + IS) = 18 -° 



1 + 3.858 x 10° pie 



5 n 2 2i<f> L 



-1 



(26) 



The corresponding result for the oscillation parameter is 



= 26.3 



1 + 3.858 x 10° pie 



(27) 
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FIG. 1: Three-dimensional plot of x s (a) and sin20 s (b) versus pi and 4>L with a Z'- mediated 
FCNC for left-handed b and s quarks. The color shadings in both plots have no specific physical 
meaning. 

With couplings of only one chirality, the physical observables AM S , x s , and sin20 s are 
periodic functions of the new weak phase <pL with a period of 180°. 

Fig. [T] (a) shows the effects of including a Z' with left-handed coupling. We see that if 
pL is small, x s is dominated by the SM contribution and has a value ~ 26. For <pL around 
90° and pL between 0.001 and 0.002, the Z' contribution tends to cancel that of the SM and 
reduces x s to be smaller than the SM value of 26.3. In Eq. (j2ZJ) and Fig.[T](a), we see that the 
Z' has a comparable contribution to the SM if p L > 0.002, independent of the actual value of 

nn 

4>l. The planned resolution of Fermilab Run II and LHCb are both about x s ^ 75 [15, 16J. 
Thus, a pl greater than about 0.003 will result in an x s beyond the planned sensitivity. If 
x s is measured to fall within a range, one can read from the plot what the allowed region is 
for the chiral Z'-model parameters. The same discussion can easily be applied to a Z' model 
with only right-handed couplings. Fig. ^ (b) shows sin20 s as a function of pl and <pL- As 
Pl increases, sin 20,, goes through more oscillations when 0^ varies from to n. 

In Fig. |21 we show the overlayed plot of the contours of fixed x s and those of fixed sin20 s . 
The shaded region in the center shows the experimentally excluded points in the <Pl~Pl plane 



9 




FIG. 2: Contour plot of x s and sm2(p s with a Z'-mediated FCNC for left-handed b and s quarks. 
The shaded region is for x s < 20.6, which is excluded at 95% CL by experiments. The hatched 
region corresponds to la ranges around the SM value x s = 26.3 db 5.5 (black curve). The solid 
curves open to the left are contours for x s = 50 (red), 70 (green) and 90 (blue) from left to right. 
The curves open to the right are contours for sin20 s = 0.5 (dotted), —0.5 (solid) and the SM value 
-0.07 ±0.01 (dashed). 

that induce x s values smaller than 20.6. The hatched area corresponds to the parameter 
space points that produce x s values falling within the la range of the SM value of 26.3. 
Contours for higher values of x s are also shown. The SM predicted sin20<, ~ —0.07 ± 0.01 
would appear as narrow bands around the sin20 s = —0.07 curves. Note that even if the 
x s measurement turns out to be consistent with the SM expectation, it is still possible that 
the new physics contributions, such as the Z' model considered here, can alter the sin20 s 
value significantly. It is therefore important to have a clean determination of both quantities 
simultaneously. Once x s and sin20 s are extracted from B s decays, one can determine up 
to a two-fold ambiguity and 4>l up to a four-fold ambiguity, except for the special case with 
sin2</> s ~ 0. 

Ar s can be determined with a high sensitivity by measuring the lifetime difference be- 
tween decays into CP-specific states and into flavor-specific states. Using the J/ip4> mode, 
simulations determine [16j that the LHC can measure the ratio AT s /r s with a relative error 
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< 10% for an actual value around —0.15. Tevatron simulations show that Ar s /r s can be 
measured with a statistical error of ~ 0.02. For a sufficiently large pL, the cos# factor in 
Eq. (|3J) increases from — 1 at <$l = 0° (mod 180°) to the maximum 1 at cj>L = 90° (mod 180°). 
We are left with the phase ranges 0° < <f) L < 30°, 60° < <p L < 120°, and 150° < L < 180° 
(mod 180°) where a 3a determination of Ar s can be made. 

Once the right-handed Z' couplings are introduced, we also have to include the new 
\AB\ = 2 operators Of R , 0\ R ^ and RR defined in Eq. (JHJ) into the effective Hamiltonian 
that contributes to the B s -B s mixing. The matrix element of O is the same as that of 
LL , while those of Of R and Of^ are 

\ 2 



[B s \0[ R \B s ) 



4 
3 



mU'LB^ R (m b ) 



m b (m b ) + m s (m b ) 



[B S \0% R \B S } = 2 



m B f B B L 2 R {m b 



(28) 
(29) 



m b (m b ) + m s (m b ) 

For the Z' coupling to right-handed currents, we define new parameters pr and the 
associated weak phase 4>r: 



pR,e 



■1<PR 



92M Z R 



(30) 



g x M z > 

At the Mw scale, we have additional contributions to the effective Hamiltonian due to the 
right-handed currents, similar to Eq. (|17j) . The terms due to the left-right mixing are 



WS ^ ^PL PR e^ +M (0^ R + RL , Oi 



(31) 



In the RGE running, the Wilson coefficient of 0^ R mixes with that of O^', the relevant 



)LR. 



anomalous dimension matrices are |4 



/ 6 



7 



(o) 



— 12 



v 

/ 137 



-6 _/V + — 

UiVc ^ N c 



7 



(1) 



^ 2AL 2 



22 

3N C 



f 



\ 



71 



- + — 

4 T iV c 



203 



200 A 7" I 

3 iVc i» 

■2 i 479 I 15 
6 "t" 2N; 



N? + ^ + & + fN c f - 



44 f 

3 " 

10 

3 



\ 



22 f 
3N C J ) 



(32) 



(33) 



"2/ "6 

where N c is the number of colors and / is the number of active quarks. At the scale of the 
B meson masses, the value of / is 5. 
We take m b {m b ) = 4.4 GeV, m s {m b 



0.2 GeV, and = 225 MeV. Following 

Mb 



ll 



Eqs. (I23M24J) . we find the effective Hamiltonian terms for the operators Off at nib to be 



wS => + of , or + Of) 




(34) 



Note that there is no contribution of the operator at the M\y scale. It is induced 
through the operator mixing in RGE running and actually has an important effect at the 
mj scale, as one can see from its Wilson coefficient in Eq. (|34l) . 

In the numerical analysis, we use the central values of B[ R (nib) = 1.753 ± 0.021 and 
B% R (nit,) = 1-162 ± 0.007 given in Ref. ^jj with the decay constant fs s the same as before. 
The predicted mass difference with all the Z' contributions included is then 



AM, = 18.0 



1 + 3.858 x 10 5 (p 2 L e 2i(t)L + p 2 R e 2i ^ R ) - 2.003 x 10 6 p L p R e ii4 ' L+4 ' R) \ ps" 1 . (35) 



The overall contribution to x, from the SM and Z' is 



x, = 26.3 



1 + 3.858 x 10 5 (p 2 L e 2i<t)L + p 2 R e 2i ^ R ) - 2.003 x lOVp^^+^J . (36) 



To illustrate the interference among different contributions, we set Pl = Pr = 0.001 and 
plot x s and sin20 s versus the weak phases 4>l and <p R in Fig. |3] (a) and (b), respectively. 

First, we note that after the RGE running the operators Of and 0% interfere con- 
structively. This can be seen from the relative minus sign between the Wilson coefficients in 
Eq. (I34J) and a corresponding relative minus sign in the hadronic matrix elements given in 
Eqs. (|28jl and (|29jl . Because of the constructive interference and the fact that the bag pa- 
rameters B^ R and are twice as large as B LL , the LR and RL operators together become 
the dominant contributions. The interference of all the terms makes x s reach its maximum 
when one of the weak phases is 180° and the other is 0° (mod 360°). If pi and p R are both 
much smaller than 10~ 3 , the variation in x s in the 'Pl-'Pr space will be indistinguishable 
from the predicted SM range. Compared to Fig. [T] (a) for Z' with only LL couplings, Fig. El 
(a) shows that even for large values of pi and p R , x s can be smaller than 20.6 due to the 
interference among all the terms in Eq. ()3fij) . The current x s > 20.6 bound excludes the 
regions with (p L + (p R ~ 0° (mod 360°). Because of the assumed equal values of pi and p R , 
the parameter space points with the same sin 20 output lie along directions that are roughly 
parallel to the <Pl + <Pr = 360° line. For the more general cases of different pi and p R values, 
the crests and troughs in Fig. El (b) are no longer parallel to the 4>l + <Pr = 360° line. 

12 




FIG. 3: x s (a) and sin20 s (b) as functions of 4>l and 4>r for pl = Pr = 0.001. The color shadings 
have no specific physical meaning. 

V. CONCLUSIONS 

In this paper we discuss the effects of a Z' gauge boson with FCNC couplings to quarks 
on the B s -B s mixing. We show how the mass difference and CP asymmetry are modified 
by the left-handed and right-handed b-s-Z' couplings that may involve new weak phases 4>l 
and 4>r. In the particular case of a left-chiral (right-chiral) Z' model, one can combine the 
measurements of AM S (or x s ) and sin20 s to determine the coupling strength (p^) and 
the weak phase (pi {4>r) up to discrete ambiguities. Once comparable left- and right-chiral 
couplings are considered at the same time, we find the left-right interference terms dominate 
over the purely left- or right-handed terms, partly because of the renormalization running 
effects and partly because of the larger bag parameters. 
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